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ABSTRACT 

We present new solar models in which we use the latest low abundances and we 
further include the effects of rotation, magnetic fields and extra-mixing processes. 
We assume that the extra-element mixing can be treated as a diffusion process, 
with the diffusion coefficient depending mainly on the solar internal configuration 
of rotation and magnetic fields. We find that such models can well reproduce the 
observed solar rotation profile in the radiative region. Furthermore the proposed 
models can match the seismic constraints better than the standard solar models, 
also when these include the latest abundances, but neglect the effects of rotation 
and magnetic fields. 

Subject headings: Sun: abundances — Sun: oscillations — Sun: interior 



1. Introduction 



The standard solar models with the latest input physics are well known to yield the solar 
structure to an amazing degree of precision, and agree with the helioseismic inversions (see 
e.g., Christensen-Dalsgaard et al. f 996; Bahcall et al. 2001). Those models use the old abun- 
dance values (Grevesse & Savual 1998, hereafter GS98). However, the standard solar models 
with the new solar mixture AGS05 (Asplund et al. 2005, hereafter AGS05) disagree with 
helioseismic constraints, e.g., the position of the base of convec tion zone (CZ) is too shallow 



and t he surface helium abundance is lower than in the Sun (jChristensen-Dalsgaard et al. 





Basu 


1998 





19911 ; iBasul Il998l : iBasu fc Antial 120041 ) . Larger discrepancie s are in the sound-speed and 



density profiles be twee n the Sun and the mo dels with low Z ( jBahcall fc Pinsonneaultl 12004 : 



Guzik et al.ll2005l ). See IBasu &: Antial ( 120081 ) for a detailed review paper. 
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A number of investigations have attempted to explain this discrepancy as a matter of im- 



creases, convective overshooting, low-Z accretion 



Christensen-Dalsgaa rd 2009; ISerenelli et al.ll2009 



Bahcall et al. 


2005 


. 2006: 


Yane & Bi 


2007: 


Guzik & Mussack I 


201ol;lTurck-Chieze et al. 



2010h and revised solar composition (Asplund et al. 2009, hereafter AGSS09). The result 



of all of these analyses show that it is difficult reproduce the helioseismic constraints with 
the standard solar model also when it includes the latest abundances. This failure could be 
due to the fact that the standard solar model neglects rotation, magnetic fields and some 
extra-mixing processes. In this letter, we show that these effects can reduce the discrepancies. 



2. Beyond the standard solar model 



Helioseismology has revealed that the Sun is rotating dif ferentially at the s urface, 
slowly in the core and almost uniformly in the radiative region ( jChaplin et al.lll999l ). The 
purely rotat i on- induced mixin g has b een considered in modelin g rotating stars, as giy en in 
Zahnl Jl992h. iMaeder fc Zahnl Jl998h . iMaeder fc Mevnetl (boooh . IPalacios et all koO$ ) and 
Mathis fc ZahiT (120041 ) . However, these models appear insufficient to reproduce the helio- 
seismically inferred internal solar rotation profile. This suggests that other effects should be 
considered in extracting angular momentum from the central core of the Sun. 

Recently, two main mechanisms have been proposed to explain the solar flat rotation pro- 
file, namely internal gravity waves (e.g. Charbonnel & Talon 2005) and magnetic fields (e.g. 
Eggenberger et al. 2005). Here, we mainly describe the efficiency of the extra -mixing caused 
by rotation and m agnetic fields, as prescribed by the Tayler-Spruit dynamo (IPitts fc Tayler 



1986; 



Spruitll2002h. The theoretica l formulation of th is dynamo is still a matter of debate 



Denissenkov &: Pinsonneaultl 120071 : IZahn et al.l 120071 ). however, lEggenberger et al.l (120051 ) 
found that the model with the Tayler-Spruit dynamo-type field successfully reproduces the 
observed solar rotation profile. Therefore, it is particularly interesting to investigate the 
effects of rotation and magnetic fields on the solar interior and global parameters. 

We present a simple scheme for dealing with angular momentum transport and element 
mixing in the solar interi or. It is based on the stellar structure equ ations which include rota- 
tion and m agnetic fields (IPinsonneault et al.l Il989t iLi et al. 1 120031 ) . The detailed derivation 
is given in lYang fc Bil (120061 120081 ). This formulation allows us to estimate the effects of 



rotation and magnetic fields on the Sun properties. The angular momentum transport and 
elements mixing can be described with two diffusion equations as follows: 



pr 



~dt 



hi 



pr 4 (D rot + D r 



dr 



(1) 
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+ 



dXi 
dt 



(2) 



where the adjustable parameters fa and fc are introduced to represent some inherent un- 
certainties in the diffusion equations. The second and third terms on the right- hand- side of 
Equation (2) are the nuclear and gravitational settling terms, respectively. In our model, 
the diffusion coefficien t D rot is associated with the rotational instability as described by 



Chaboyer et al.l (119951 ). In the case of a Tayler-Spruit dynamo- type field, the diffusion c oef- 



ficient for the angular momentum transport can be written as (IMaeder fc Meynetll2003l ): 



D m = r 2 ttq 2 
and the one for chemical element transport as: 

D' = rW 



(3) 



(4) 



Equations (3) and (4) are valid in the regime of negligible thermal diffusion, namely ^> 
Nt, where Nt {N^) represents the thermal (//-) gradients associated with buoyancy fre- 
quency. When this condition is violated, we should replace Equations (3) and (4) with 



D, 



r 2 Q — 



n 

N T 



1/2 



K 



and 



respectively, where q 



D' 



91nQ 



r 2 Vt \q\ 



iM 3/4 

Nt) 



r 2 N 7 



K 

r 2 N 7 



1/2 



(5) 



3/4 



(6) 



n and K = AacT 3 /3np 2 c p is the thermal diffusivity. 
Add i tiona lly, in order to reproduce the solar surface angular velocity, we adopt the 



Kawaler Jl988h braking law: 



dJ 
~dt 



ficK c 



R 

Rq 



1/2 



M 



-1/2 



Q 3 



(7) 



where ~ 1.13 x 10 47 g cm 2 s and fx is an adjustable parameter related to the magnitude 
of the magnetic fields. 
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Calculations and results 



Our solar models are obtained from the one-dimensional Yale Rotating Stellar Evolution 
Code (YREC; Guenther et al. 1992; Li et al. 2003; Yang & Bi 2006), by including rota- 
tion, magnetic fields and relevant extra-mixing processes. In addition, we use the following 
upda ted physical quantities: OPAL equation of state tables EOS2005 ([Rogers fc Nayfonov 
20021 ) , the opacities (GS98, AGS05 and AGSS09) supple mented by the low- temperature opac- 
ities (IFerguson et al.l 120051 ). diffusive element settling (IThoul et al.lll994j ) and the Krishna- 
Swamy Atmosphere T — r relation. 

In order to investigate the influence of rotation and magnetic fields, we constructed 
these solar models in accord with different physical processes corresponding to different 
solar compositions. In the numerical calculations, all models are calibrated from the initial 
zero-age main sequence (ZAMS) to the present solar-age models, for which the radius is 
6.9898 x 10 10 cm, the luminosity 3.8515 x 10 33 erg/g, the mass 1.9891 x 10 33 g and the 
adopted photospheric Z/X ratio. The free variables are the initial helium abundance Y, the 
initial metallicity Z and the mixing-length parameter, all of which are adjusted to match 
these observational con straints. In addition , we a ssumed that the convective region rotates 
rigidly, as proposed by iPinsonneault et al.l (119891 ) . The initial angular velocity is another 
free parameter that can be tuned so that the surface velocities of solar-age models match 
the observed values. 

Figure 1 shows the angular velocity as a function of radius r at the ages of 2.0 Gyr and 
4.57 Gyr. For the purely rotating model, the ^-gradient clearly appears in the radiative 
region. In the solar interior, the angular velocity increases with increasing age during the 
main sequence, while in the surface it is just the opposite. As a consequence of this at the 
present age the core rotation velocity is about four times as large as the surface one. On the 
other hand, the angular velocity profile for the model with magnetic fields is significantly 
different. During the main-sequence stage, the Sun is a quasi-solid body. It is interesting to 
note that at the age of 4.57 Gyr, the surface rotation velocity predicated by both models is 
approximately 2.9 x 10 -6 rad/s. However, the total angular momentum is quite different. 
For the calibrated models, the total angular momentum of the rotating model at the age 
4.57 Gyr is 8.97 x 10 48 g cm 2 s " 1 , which is about fiy e times as large as the seismic result 
(1.94 ± 0.05) x 10 48 g cm 2 s" 1 flKomm fc Howd 120031 ): while for the model with magnetic 



fields at the same age, the total angular momentum is 2.02 x 10 g cm s , which is in good 
agreement with the result obtained by helioseismology at la level. The magnetic field thus 
constitutes a more efficient process to transport angular momentum because it enhances the 
coupling between the radiative zone and convective one. 



Rotation and magnetic fields have important consequences on the chemical composition 
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profile of the outer convective envelope, as shown in Figure 2. By comparing the models with 
and without rotation and magnetic fields, we find that the extra-mixing process counteracts 
the effect of diffusive settling in the outer envelope. Hence, the model with rotation and 
magnetic fields has a smoother helium abundance profile than the standard solar model with 
the same abundance. This leads to a change in the CZ structure, as well as improvements 
in the sound speed and density profiles. 

For further investigation of the magnetic field's role, Figure 3 shows i n detail the differ - 



ences between the calculated and inferred sound-speed and density profile (IBasu et al.ll2009l ). 
The different lines refer to calibrated evolved models at the age of 4.57 Gyr, each using a 
different abundance, indicated as GS98, AGS05 and AGSS09. It is clearly visible that the 
discrepancy between seismic inferences and solar models with the new lower abundances is 
much larger than one with the old abundance. As illustrated by the different curves, the 
effects of rotation and magnetic fields on the stellar structure equations change the hydro- 
static equilibrium and thermodynamic variables on the solar interior, and therefore also have 
a significant impact on the solar models. Table 1 summarizes the main characteristics of our 
calibrated models. Interestingly, among all the models listed in the table, model AGSS09c 
shows the best agreement with the inversions. This model reproduces the sound speed and 
density profiles in within 0.5%, while for model GS98a the discrepancy is about 0.3%. Fur- 
thermore, model AGSS09c predicts the position of the CZ base at Rcz = O.721i? which 
shows a 8(T discrepancy, while AGSS09a model shows a 10a one. For the surface helium 
abundance the situation is analogous: model AGSS09c predicts Y s = 0.243 with a 1.5a dis- 
crepancy, while model AGSS09a shows the discrepancy at 3.6a level. Although the models 
including rotation and magnetic fields show some improvements with respect to the standard 
solar model, they still disagree with the seismic constraints. 



4. Conclusions 

We have investigated the effects of rotation and magnetic fields on the solar models, and 
have found that when these effects are included, alongside the new abundances, the revised 
solar models can better reproduce the helioseismic constraints. However, we see that it is 
difficult to match simultaneously the new abundances and helioseismology data for sound 
speed, density profiles, convection zone depth and surface helium abundance. Although 
the Tayler-Spruit dynamo type magnetic field still needs to be studied further, our results 
show that it does provide a possible explanation for the solar abundance problem. We have 
neglected turbulence, which may feed the differential rotation and sustain magnetic fields 
in the convection zone, and other interactions. These physical processes will be considered 
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in our future work. The results obtained in this paper are encouraging and we intend to 
apply our model to solar-type stars to get a proper interpretation of the existing helioseismic 
observations and the coming asteroseismic ones. 
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Table 1. Characteristics of the Calibrated Solar Models 



Model 


(Z/X) s 


Z 8 


Y s 


Rcz/Rq 


< 5c/ c > 


< Sp/p > 


Y c 


Z c 


Y 




OtMLT 


GS98a 


0.0229 


0.0169 


0.246 


0.715 


0.0012 


0.008 


0.644 


0.0198 


0.277 


0.0188 


2.12 


AGS05a 


0.0165 


0.0125 


0.230 


0.728 


0.0030 


0.034 


0.623 


0.0148 


0.261 


0.0140 


2.08 


AGS05b 


0.0165 


0.0124 


0.239 


0.727 


0.0028 


0.035 


0.622 


0.0146 


0.269 


0.0139 


2.04 


AGS05c 


0.0165 


0.0124 


0.237 


0.726 


0.0028 


0.033 


0.621 


0.0146 


0.260 


0.0139 


2.05 


AGSS09a 


0.0181 


0.0136 


0.236 


0.723 


0.0020 


0.022 


0.631 


0.0160 


0.267 


0.0152 


2.12 


AGSS09b 


0.0181 


0.0134 


0.245 


0.722 


0.0019 


0.023 


0.630 


0.0158 


0.268 


0.0150 


2.07 


AGSS09c 


0.0181 


0.0135 


0.243 


0.721 


0.0017 


0.021 


0.630 


0.0158 


0.266 


0.0150 


2.09 



a Solar models with diffusion. 

b Solar models with diffusion and rotation. 

c Solar models with diffusion, rotation and magnetic fields. 
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Fig. 1. — Comparison of angular velocity profiles for the two cases: (1) purely rotating 
models; and (2) models with rotation and the Tayler-Spruit dynamo-type field. The dashed 
and solid lines refer to different ages: 2.0 Gyr and 4.57 Gyr, respectively. 
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Fig. 2. — Helium abundance profiles for the calibrated solar models at the age of 4.57 Gyr, 
computed with and without rotation and magnetic fields. 
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Fig. 3. — Differences between inferred and calculated sound speeds and densities for models 
with and without rotation and magnetic fields at the age of 4.57 Gyr, corresponding to the 
GS98, AGS05 and AGSS09 abundances. Sound speed and density inversions are from Basu 
et al. (2009). 



